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Abstract 

We report theoretical results demonstrating that photoionization can be a useful tool for investi- 
gating single-electron transistors. It permits to obtain information on the quantum dot occupancy 
and the charging energy in a direct manner, and not indirectly, as done in transport experiments. It 
is worth emphasizing that in the photoionization processes considered by us, an electron absorbs a 
photon with energy of the order of the work functions (~ 1 eV) and is ejected into the vacuum. This 
phenomenon is completely different from the widely investigated photo-assisted tunneling. There, 
an electron tunnels through a Coulomb island from one electrode to another by absorbing a photon 
of much lower energy, of the order of the charging energy (typically, a few meV) . Suggestions are 
given on how to conduct experiments using photoionization alone or in combination with transport 
measurements. Monitoring zero kinetic energy (ZEKE) photoelectrons is especially recommended, 
because ZEKE-spectroscopy offers a better resolution than standard photoemission. 
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I. INTRODUCTION 



In a single-electron transistor (SET), a quantum dot (QD) is coupled via tunneling 
junctions to metallic electrodes An isolated quantum dot (QD) behaves like an 

"artificial" atom, where the electron motion is confined within a small region with size of 
the order of nanometers ^, ^J. Similar to atoms, the single electron levels of a QD are 
well separated energetically. Experiments on SET transport are often interpreted within 
the Anderson impurity model at equilibrium The QD is modeled by a single 

localized state. Its energy Ed can easily be tuned by varying the potential V g of a gate. 
This approximation is particularly justified in small QDs, where single electron levels are 
well separated energetically. This state can accommodate = 0, 1, or 2 electrons. A key 
quantity for electronic transport through the SET is just this number rid of ("valence") 
electrons. 

Unlike in extended systems, a nonvanishing charging energy U = e 2 /2C has to be paid to 
add one electron to a QD (C is the total capacitance). At sufficiently low temperatures, for 
a small QD weakly coupled to electrodes, the charging energy U represents the dominant 
energy scale; it exceeds the thermal energy and the finite width T of the dot level resulting 
from its hybridization (td) with electron states of the electrodes. 

Therefore, charge fluctuations are largely suppressed and the number of electrons of the 
dot is an integer in broad e^-ranges (plateaus): rid ~ for e d > e f + V, rid — 1 for 
Ef — U + T < Ed < £f — r, and rid — 2 for Ed < £f — U — V. The transitions between 
the states with well defined dot charge occur within narrow e^-ranges of width ~ 2T. In 
these ranges, the QD state is a combination of two nearly degenerate states of well defined 
charge (mixed valence regimes). There, the dot charge can fluctuate and, thus, electron 
transport becomes possible. This gives rise to the well known Coulomb blockade peaks, 
separated by Aej — U in the curve of the conductance G(ea) at low temperatures (but 
larger than the Kondo temperature Tk) [4]. Although charge fluctuations are suppressed 
within the plateaus, spin fluctuations through virtual intermediate state are allowed within 
the (Kondo) plateau where the QD is occupied by one electron (rid — 1)- At T < T K , they 
yield a narrow peak (Kondo resonance) in the QD density of states, which enhances the 
conductance up to the value of the ideal point contact Gq = 2e 2 /h (unitary limit) 

The experimental observation of the Coulomb blockade phenomenon and the Kondo effect 
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in electric transport of SETs is remarkable 0, [lo|. It indirectly confirms the occurrence 
of the charge plateaus in broad V^-ranges. This evidence is indirect, because what one 
directly measures in experiment is the conductance G, and not the dot charge n^. Besides, 
the charging energy U is also indirectly estimated as Asd ~ U, while the experimentally 
measured quantity is V g and not 

In the present paper, we shall propose an alternative method to investigate phenomena 
related to the charge plateaus in a SET, namely the photoionization. Photoionization is 



known as a very useful tool to study strong electron correlations [ll|, LL2J- As we shall show, 
from the investigation of the ionization one can extract more direct information on the QD 
charge than from transport experiments. 

The remaining part of the paper is organized as follows. In Section [III we describe the 
model and present the theoretical results on photoionization. Based on these theoretical 
results, in Section [Till we suggest possible experiments to employ photoionization as a tool 
for investigating SETs. Conclusions are presented in the final Section | 

II. MODEL AND THEORETICAL RESULTS 

We shall also describe the SET by the Anderson model 
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where ai >a {a\ a ) denote annihilation (creation) operators for electrons of spin a in the left 
and right leads (L,R), d a = ao, CT (d-l = Oq jCT ) destroys (creates) electrons in the QD, and 
n^a = d) a d a are electron occupancies per spin direction. The QD-electrode coupling is 
characterized by the hopping integral td, which can be experimentally controlled. In the 



experimental setup of Ref. 
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ta can be changed by varying the potential of the gates that 
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form the constrictions, and the dot energy Ed can be tuned by varying the potential of a 
"plunger" gate electrode V g 

Ed = aV g + const. (2) 

The number of electrons will be assumed equal to the number of sites (N = Mi + Mr + 1). 
In the photoionization of the QD, a photon of energy u is absorbed by an electron, which is 
ejected from the QD into vacuum, and brings the ionized system in one of its eigenstates ^k- 
The ionization process can be characterized by an energy threshold Uk=(^f k\H\^f k)-(^\H\^) 
and a spectroscopic factor f^ a 

fk,a= \{* k \dM\ 2 - (3) 

Here $ stands for the neutral ground state (case T = 0). The spectroscopic factor is 
directly related to the weight of a line in the ionization spectrum, more precisely to the 



ll| . From Eq. (j3J) one can deduce the following sum 



partial-channel ionization cross section 
rule 

Y, h,a = ($14^1$) = n d,a- (4) 
k 

Eq. (j4j) represents an important result. It permits to directly relate the integrated weight of 
the ionization spectrum to the number of electrons on the QD in the neutral ground state. 
From a technical side, Eq. (T4J) turns out to be very useful to test the numerical results. It 
allows to check whether the Lanczos algorithm targets all the relevant ionized states ^k- 

To compute the ionization spectrum, we shall employ the method of exact numerical 
(Lanczos) diagonalization. Although this method can only be applied to small electrodes, 
the results are relevant provided that the number of sites of the electrodes are properly chosen 
to mimic a metallic behavior jl3l |. The single-particle energies of an isolated electrode lie 
symmetrically around ep within the range {ep — 2t, Ep + 2t). To ensure that in the ground 
state of the isolated electrodes the Fermi level is occupied by one electron, one should 
consider that each electrode consists of an odd number (M L p) of sites. To demonstrate the 
reduced role of the finite-size effects, we have also considered the QD embedded in a ring 
(a-M L ,a = a M R ,a)- In this case, for a similar reason, the ring without the QD should consist 
of an odd number of sites (odd Mp + Mr). In all numerical results presented here, we set 
t — 1. For the numerical calculations, the value of the Fermi energy in electrodes Ep is not 
important: it only fixes the energy zero. Therefore, we do not fix the value of Ep and present 
in all the figures results for energies relative to e_p. 
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Our numerical results for rid of Fig. [T] demonstrate the formation of broad, well defined 
charge plateaus (rid ~ 0,1,2), separated by narrow mixed valence regions, as discussed 
above. The fact that the n^-curve (as well as those for the bright diabatic states, see below) 
for various electrode sizes and boundary conditions are almost identical within the drawing 
accuracy in Fig. [T] demonstrates that finite-size effects are not important and the study of 
small systems is physically relevant. Our extensive numerical results demonstrate that what 
is essential for the weak finite-size effects is also essential for the occurrence of the charge 
plateaus, namely a width parameter T oc t\jt smaller than the charging energy U. On this 
basis, we argue that the charge and the ionization spectra of a QD connected to semi-infinite 
leads can be accurately deduced by computations for such small "metallic" electrodes. 

Although it does not represent a central issue from the present perspective, the weakness 
of finite-size effects is also relevant for the Kondo effect, which we only mention here in 
passing. As is well known, essential for the occurrence of this effect is the Kondo screening 
cloud. The latter extends over a number of sites £x ~ t/Tjc, which rapidly grows (exponen- 
tially for large U) far beyond electrode sizes treatable by exact numerical diagonalization. 
Therefore, it is impossible to obtain the conductance within standard transport approaches 
by employing short electrodes. For semi-infinite electrodes, the zero-bias conductance G can 
be computed via the Friedel-Langreth sum rule Q, 16] 



G = G sm 2 (7rnf/2) (5) 

once the dot occupancy nf is known. Eq. dSJ) was deduced for a QD attached to semi- 
infinite leads with a continuum density of states, which is the real nanosystem of interest. 
The pleasant thing is that the zero-bias conductance G for this case can be solely determined 
from n^°, and a reliable estimated of the latter suffices for this purpose. Fortunately, to this 
we need not to carry out calculations for semi-infinite or very large electrodes: in view of 
the weak finite-size effect mentioned above, rid computed for short electrodes does represent 
an accurate estimate, nf ~ n d . Indeed, the curve for G(ed) obtained by inserting the 
calculated rid calculated for short electrodes instead od n^°, also presented in Fig. [lj has 
a similar appearance to that obtained by means of the numerical renormalization group 
171 . Il8j ] for semi-infinite electrodes and nicely reveals the occurrence of the Kondo plateau. 

Let us now switch to the results for SET photoionization. To avoid misunderstandings, 
we emphasize that the photoionization considered in this paper is the qualitatively different 
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FIG. 1: (Color online) Results for the dot occupancy rid and normalized conductance G/Gq for 
U = 8, td = 0.2, t = 1, and Ml = Mr = 5 (N = 11). These curves can be hardly distinguished 
within the drawing accuracy from those for open chains with (Ml, Mr)=(3, 3), (3,5), and (3,7), 
as well as for rings with N = 6, 10. 
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from the widely-studied photo-assisted tunneling (see, for example, Ref. 
therein). The process where a photon helps an electron to tunnel from one electrode to 
another through a Coulomb island might also be viewed as an "ionization" phenomenon. 
Nevertheless, the scale of the corresponding "ionization" energies is determined by the charg- 
ing energy U. For semiconducting QDs, it is of the order of ~ 1 meV, that is, much smaller 
than ionization energies of the order of the work functions (~ 1 eV) considered by us. 

Typical exact numerical results for the ionization energies and spectroscopic factors of 
all the significant ionization processes are depicted in Fig. [2j To understand these results it 
is helpful to analyze first the limit — >• 0. 

In this limit, two ionization processes of the QD are possible. For Ed < Ep — U, the 
dot is doubly occupied (n<j = 2), and the relevant ionization process consists of removing 
an electron from the upper Hubbard "band" . The corresponding ionization energy is = 
—Ed — U. For Bp — U < Ed < £f, the dot is occupied by a single electron in the lower 
Hubbard "band" (rid — 1)- To remove it, an ionization energy uf = is needed. The 
corresponding spectroscopic factors are /° = 1 and /° = 1/2. 

Basically, the exact results presented in Figs. [2^i and [2b differ from those of this limiting 
case in two ways. First, in the mixed-valence regimes Ed ~ ef and Ed ~ Ef — U the finite 
coupling causes a smearing effect similar to that already observed on the curve nd(£d) in 
Fig. [TJ Second, the exact curves exhibit numerous avoided crossings. In the present case, 
avoided crossings occur around the points where elementary ionization processes become 
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quasi-resonant. All nearly linear pieces of the curves in Fig. [2^ are reminiscent of elementary 
ionization processes in the noninteracting case (td —>■ 0). For instance, one-hole processes 
in electrodes are represented by horizontal lines. Two-hole-one-particle (2h-lp) processes 
where the creation of one hole in the upper or lower Hubbard band is accompanied by the 
excitation of a particle-hole pair in electrodes are represented by lines parallel to the u®- 
or oof- lines, respectively. To mention only one more example, 2h-lp processes involving the 
creation of two holes in the upper and lower Hubbard bands and of one electron in electrodes 
at or above the Fermi level are represented by the straight line uf, = —2e d — U + ep or 
parallels to it. At the intersection points of such lines avoided crossings arise rather than 
true intersections, the latter occurring only for t d — > 0. As discussed recently [20|, avoided 
crossings are very frequent in tunable QD systems. 

Numerous avoided crossings are visible in Fig. [2^, but for practical purposes they are less 
important, and therefore only a few are shown in Fig. [2b. Parts thereof are not important 
because the spectroscopic factors are too small. But even if the signals are sufficiently 
intense, avoided crossings are less important: in the narrow regions around the critical 
points where the intensities are significant, the energy differences are usually smaller than 
the experimental resolution, and only the summed intensities of the participating states can 
be measured. What one measures is a smoothly varying spectroscopic factor, attributed to 



a bright "diabatic" state 



20] . The constant spacings between different families of parallel 



lines in Fig. [2k can be understood by noting that the single-particle energies of the isolated 
electrodes considered there are £p, Sf ± t, and Sp ± t\f?>. For long electrodes, bundles of 
dense parallel lines within a width At will appear, so of practical interest are just such bright 
diabatic states passing through a multitude of avoided crossings. 

Therefore, besides the results for individual eigenstates, we also show in Figs. [2^ and 
[2b the curves Q u j and f u j for the two bright diabatic states that are significant. Deep 
inside the plateaus, the ionization signals corresponding to these two diabatic states behave 
similarly to the case td — > 0, but substantial deviations are visible in the mixed valence 
regimes. Noteworthy is not only that the spectroscopic factors penetrate the neighboring 
charge plateaus, where they gradually decrease, similar to n d (Fig. [I]). Even more interesting 
is the fact that all these penetrations are accompanied by substantial changes in the slope 
of the fl u /-curves. 



7 











(a 


) 


















1 X. \ ' "X 2 



-12 -8 -4 e,-e f 4 









(b) 






2 






■ 






































1 1 









-12 -8 -4 e.-e_ 4 

a f 



FIG. 2: (Color online) Results for ionization energies (a) and spectroscopic factors (b). The solid 
lines 1, 2, 3, and 4 are for the eigenstates with substantial spectroscopic factors. The results for 
the two relevant diabatic bright states (u, I) are depicted by symbols. The black points in (a), 
which correspond to the ionized eigenstates with small spectroscopic factors (10~ 5 ^ / i$ 10~ 2 ), 
are not quantitatively significant for the ionization spectrum, but are showed to better visualize the 
occurrence of the avoided crossings. Notice that the ionization energies are measured relative to 
the electrode work function (— ep), a quantity much larger than the charging energy U. Parameter 
values as in Fig. [TJ See the main text for details. 

III. PROPOSED EXPERIMENTS 

In view of the present theoretical results, we propose to conduct the following experiments 
using an intense incoming flux of monochromatic photons, which is well focused on the QD 
of a SET. To avoid confusions we emphasize again that the photons should have energies 
of the order of the work functions, significantly higher than of the microwave and rf-fields 
considered in the previous studies on SETs and the Anderson model. 

(i) The first type of experiments is a standard ionization study. To measure the absorption 
coefficient would be desirable but presumably a too hard task. Probably it would be easier 
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to measure the zero kinetic energy (ZEKE) photoelectrons emitted from the dot at the ion- 
ization threshold, as in threshold ionization experiments in molecular physics. Just because 
the QD is very small, not only photoelectrons from the QD, but also from electrodes will be 
inherently ejected into vacuum. Therefore, it is important but fortunately easy to experi- 
mentally distinguish between the ionization of the QD and that of the electrodes: by varying 
the gate potential V g , the ionization energies u>i change in the former case [ui ~ — Sd(V g ) 
and u)i ~ —EdiVg) — U] but remain constant in the latter {oj-i = — ep). Let us analyze the 
change in ionization by gradually increasing the dot energy e d starting from a sufficiently 
low value, by making V g more and more negative. In the first stage, an ionization signal 
will be observed, which is characterized by a nearly constant intensity I u {I u oc w u ~ 1). Its 
intensity will drop to zero within a narrow range SV g i ~ 2AV (AV = —T/a > 0) centered 
on V 9 = V g \. A little before this signal disappears, another ionization signal will appear at 
V g ~ V g i + AV, whose intensity rapidly rises to a value Ii « I u /2 (Ii oc wi ~ 0.5) beyond 
Vg ~ Vgi — AV. Further on, this intensity remains nearly constant up to a point where it 
rapidly falls down to zero within a range 5V g 2 ~ 2AV centered onV g = V g 2- Except for the 
(mixed-valence) ranges 5V g i and SV g 2, the corresponding ionization energies, Q u (V g ) and 
fl[{V g ), are straight lines of the same slope. This slope is half the slope of the f2/(V^,)-curve 
in the range V g i < V g < V g \ + AV. At the other end, V g 2 — AV < V g < V g 2, this curve 
tends to saturate, VLi(V g ) — > —ep. The same tendency of the other curve (fl u (V g ) — > —ep) 
can be seen in the opposite mixed-valence range V g \ — AV < V g < V g \. By extrapolat- 
ing the linear parallel portions of the f2 Uj ;-curves (Q u ,i —* &>u,i), one can directly determine 
the charging energy from the difference of their ordinates taken at the same value of V g , 
U = Q u {V g ) — VLiiVg) [2l|. Noteworthy, for this, the values of V g or e d themselves are not 
needed. This represents an important difference from the determination of U in transport 
experiments. In the latter, what one can directly measure are the values V g \ and V g 2, from 
the positions where the Coulomb blockade peaks occur at temperatures T > Tk or from 
the extension of the Kondo plateau observable for T < Tk- The determination of U from 
the difference V g i — V g 2 alone is impossible, and to this aim the relationship between V g 
and Ed, Eq. (j2J), is necessary. There is still another important difference between ionization 
and transport measurements. As visible in Fig. [2^, fl u (V g i) ~ fi^V^) — Sp. The Fermi 
energy can thus be estimated. This is also important because this quantity can alternatively 
be deduced from the ionization energy — ep of the electrodes (electrode work function). To 
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conclude, the ^-dependent ionization spectrum allows one to give direct evidence on the 
formation of the charge plateaus and to determine the quantities U, T, and Sp. The depen- 
dence SdiVg) can also be obtained, which is also important, as it allows to quantitatively 
investigate the capacitances of nanosystems. 

(ii) A (ZEKE) ionization study might be difficult, but a mixed ionization-transport ex- 
periment as proposed below is presumably easier while still relevant. In a usual transport 
experiment, one can measure the conductance as a function of the gate potential, G(V g ), 
and determine V gl and V g2 , say, from the Kondo plateau (T < T K ). At a fixed gate potential 
Vg2 < Vg < Vgi one can now shine the QD with radiation whose frequency u is varied. At a 
certain value of the latter, uo = uj c (V 9 ), the incoming photon will ionize the QD. Because the 
single (unpaired) electron of the QD will be ejected, the Kondo effect will be suppressed, 
and this will be evidenced by the drop in the conductance (G ~ 0). From the above consid- 
erations one expects a linear V^-dependence of cu c (V g ) = Q[(V g ). Again, the charging energy 
can be directly determined as U = £li(V g i) — Q;(V^ 2 )- The dependence of Ed on V g is not 
needed, rather it can be deduced by this method, as well as the Fermi energy e F = —Vti(V g 2). 

(iii) A similar mixed experiment can be imagined at V g < V g 2, where the nonionized 
QD is doubly occupied. At the resonant photon frequency u = Q u (V g ) one electron will be 
removed. For a sufficiently intense irradiation, one can reach a stationary state where the 
ionized QD will be occupied by an unpaired electron. Because the latter is a prerequisite for 
the occurrence of the Kondo effect, it is tempting to speculate on a possible photoionization- 
induced nonequilibrium Kondo effect; thence, an effect of photoionization opposite to that 
of suppressing the equilibrium Kondo effect discussed above. 

IV. CONCLUSION 

In this paper, we have showed that photoionization is a valuable method, which can be 
used along with transport measurements, to study SETs. To our knowledge, conductance 
measurements represent the only way to determine the charging energy U, which thus resem- 
bles very much a fit parameter. Therefore, the comparison with the result of a completely 
different method is quite meaningful, even more if U is obtained in a direct manner, and 
not via a supplementary relation between Ed and V g . For specific purposes, photoionization 
turns out to provide even richer information than transport studies. 
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The above considerations refer to a QD with a single level. However, from the physical 
analysis it is clear that the results are also relevant for more general situations. 

Another important result is that QD-ionization can be accurately described by consider- 
ing short electrodes (weak finite-size effects), and only very few ionized states are relevant. 
Both aspects are very important, since they urge to carry out realistic ab initio calculations 
on photoionization, particularly for single-molecule based SETs. 
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